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Abstract 

We consider a one-dimensional mesoscopic Hubbard ring with and without 

disorder and compute charge and spin stiffness as a measure of the permanent 

currents. For finite disorder we identify critical disorder strength beyond 

which the charge currents in a system with repulsive interactions are larger 

than those for a free system. The spin currents in the disordered repulsive 

Hubbard model are enhanced only for small U , where the magnetic state of 

the system corresponds to a charge density wave pinned to the impurities. 

For large U , the state of the system corresponds to localized isolated spins 

and the spin currents are found to be suppressed. For the attractive Hubbard 

model we find that the charge currents are always suppressed compared to 

the free system at all length scales. 
72.15Rn, 73.20Mf 
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I. INTRODUCTION 



Over the last few years various experiments have measured the magnetic response of 
small, moderately disordered ensembles of quasi one-dimensional (ID) rings.0^ These ex- 
periments confirm the existence of a persistent current which had been predicted already 
much earlier.0 However, the experimental value for the persistent current is two to three 
orders of magnitude larger than the theoretical predictions based on calculations in a disor- 
dered but non interacting electron gas.i It is thus commonly believed that the discrepancy 
could probably be resolved by accurately including the electron-electron interactions in these 
calculations. 

Studying the full interacting electron problem with disorder is in general difficult. How- 
ever, in one dimension, powerful analytical tools, such as Bethe Ansatzifl and Bosonizationi'i 
give us a handle to treat at least the on-site part of the interaction exactly and the ef- 
fect of disorder can then be studied within a perturbative renormalization group (RG) 
approach. Using these methods, Giamarchi and Shastrylll have shown that a repulsive inter- 
action enhances the value of the persistent current in a mesoscopic Hubbard ring. This result 
has been confirmed by recent independent numerical^ and analytical studies treating the 
Hubbard interaction in first order perturbation theory. 

In this paper we study the persistent currents in the ID mesoscopic Hubbard ring both 



with and without disorder along the lines of Ref. |I0|. Due to the absence of Galilean 
invariance,lli we note that the persistent currents exhibit a strong reduction as a function 
of both filling and interaction strength already for the clean ring.0 Thus it is a priori not 
clear that this initial reduction can be compensated completely by disorder as to give a net 
enhancement of the current. Therefore, we use the exact bosonization parameters, obtained 
from the Bethe Ansatz solution for the clean ring, as initial values for the solution of the 
RG equations. We remark that Ref. ^ uses only perturbative starting values for the RG 
equations and normalizes the current w.r.t. non-interacting values. Our approach allows us 
to show that for arbitrarily weak disorder there is no enhancement in the currents. Fur- 
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thermore, we can identify critical disorder strengths beyond which the persistent current 
in the presence of both interactions and disorder is indeed larger than the current for the 
disordered but non interacting system. 

We also study the effect of disorder on the spin currents in the Hubbard ring. We 
emphasize that even in the presence of Galilean invariance the spin currents can renormalize 
non-trivially under interactions. Physically a non-vanishing spin current implies that the 
system has long ranged spin cor relat ions, and it is thus interesting to see the effect of 
disorder on the spin currents. 

The plan of the paper is as follows: In section |I| we define charge and spin stiffness for 
a ID mesoscopic system and relate them to the corresponding currents. We briefly discuss 
the range of validity of these relations. In section |T| we first review the Bethe Ansatz (BA) 
solution for both the repulsive and the attractive clean Hubbard model. We then go on to 
compute the stiffnesses and their conjugate compressibihties by iterating the BA equations 
for a mesoscopic system. Next, we recall the bosonized description of the Hubbard model 
and show how to relate stiffness and compressibility to the new parameters of the bosonized 
Hamiltonian. In section |rV|, we then include the effects of disorder into the bosonized model 
by a RG calculation. These RG equations are then integrated numerically and we can 
identify a crossover from a phase with enhanced spin currents to a phase in which spin 
currents are suppressed as the on-site interaction is increased relative to the initial disorder 
strength. We discuss these results in section 0. 

II. PERSISTENT CURRENTS AND STIFFNESS OF A ID MESOSCOPIC 

SYSTEM 

For a finite system on a ring of length L, the response of the ground state energy to 
a finite Aharanov-Bohm (AB) flux is a measure of the persistent current at T = 0. The 
current is given as J = L (9i?o($)/c^$|$=o' where Eo{^) is the ground state energy of the full 
interacting system in the presence of the AB flux $. The energy shift of the ground state 
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can be written as Eo($) - £^o(0) = D^'^ / L + 0($4), where D = {L/2)d'^Eo{<^)/d<^^\^^^^ is 
called the stiffness constant. D provides an operational definition for the persistent current 
for small values of the flux, given as J = 2D^. Higher order terms are important when the 
energy shift is comparable to the mean energy level spacing in the spectrum of the many- 
body system. In the case of a finite ID metallic system, the gaps are 0(1/ L) and so non- 
quadratic corrections occur when $ is 0(1). Level crossings would occur and perturbation 
theory would break down for $ of order vr. 

In a system with up and down spins, the AB flux for each species can be treated as 
independent parameters $| and $j. With this freedom, there are two stiffnesses that can 
be calculated and hence two currents. In the case when $| = $| = $c, both species are 
coupled to the same flux and the shift in the energy gives the charge stiffness Dc- In a 
Galilean invariant system this flux couples only to the center of mass coordinate and hence 
the persistent charge current will not depend on the interaction between the particles.lli 
However, for lattice models such as the Hubbard model, the center of mass momentum is 
conserved only up to the reciprocal lattice vector and we can observe a non-trivial dependence 
of the charge current on the interaction. In the sector of zero magnetization {Sz = 0), the 
second choice is to take $| = — $1 = $s/2. In this case the two species are driven in opposite 
directions through each other which leads to a non-trivial dependence of the resulting spin 
current on the interactions between opposite spins even for a Galilean invariant system. The 
energy shift in the presence of $s defines the spin stiffness D^. Charge and spin currents are 
then given as Jr = 2Dj.^r and 

Ld^Eoi^r 



Dr 



(1) 



We note that Dc is related to the d.c. part of the conductivity cr(co') = D£,rude^{.^) + c^reg by 



III. THE ID HUBBARD MODEL WITHOUT DISORDER 
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A. The repulsive interaction regime 



1. Charge and spin stiffness and the conjugate compressibilities 

The repulsive Hubbard model on a ring of size L threaded by a spin dependent flux $0- 
is described by the Hamiltonian 

H = -tj2 (e^*^^''4+i,.Q,^ + h.c) + uj^n^^^n,^^, U > 0, (2) 

i=l i=l 

where c- g. and q^o- create and annihilate fermions on site i with spin a and ni^„ = c\^„Ci^a is 
the density operator. We define our energy scales by setting t = 1 in the following. Note 
that this is in fact equivalent to using the Fermi velocity of the non interacting system Vf 
as defining energy scales, as f / = 2tsin[7™/2], where n is the density. 

The addition of a flux $0- is compatible with integrability and the Bethe Ansatz equations 
for a chain with total number of particles N = + Ni and M = number of down spin 
particles are given aS 

M 

Lkn = 2nln + $1 + 2 ^ arctan [4(Aj - sin kn)/U], (3a) 

AT M 

2 arctan [4(Aj - sin /U] = 2TiJj + $| - $1 + 2 ^ arctan [2(Aj - A^) /U]. (3b) 

n=l j^j 

The quantum numbers /„ and Jj label different states and statistics. For the fermionic 
ground state they are taken to be A/„ = M/2 mod 1 and AJj = (A^ — M + l)/2 mod 1. 
The energy of the system in a state corresponding to a solution of Eq. (|^) is equal to 
Eo{^^) = -2Ecosfc„. 

In order to study the currents in the ring one should in principle distinguish between 
cases of odd or even numbers of particles: If both N-^ and Ni are odd, the number of right 
and left moving particles for each of the species are equal, the ground state is non-degenerate 
and the energy is minimum in zero external flux. For any other combination, i.e. when either 

or Ni or both are even, the number of right and left moving particles are not the same 
due to the /c = state and hence the energy is not necessarily a minimum at zero charge or 



spin flux but rather at flux values of ±7r. For the sake of simphcity, we shall only consider the 
case of zero magnetization with N even and M = N/2 odd in our numerical computations. 

In Figs. and ^ we show Dc and Dg as obtained from an iteration of the BA equations 
for a ring of size L = 100 and various flllings n = N/L and values of the interaction strength 
U. Both figures show that there exists a reduction in both stiffnesses with increasing U for 
a given filling. The magnitude of the reduction for Dc is considerably enhanced as the filling 
is increased towards the half-filled case. 

The other physical quantities that will be needed in the next section to compute the 
parameters arising in the bosonized Hubbard Hamiltonian are the charge compressibility and 
the spin susceptibility. The usual thermodynamic expression for the charge compressibility 
is given as Xc^ ~ {''^'^ / L)d'^EQ/dn'^\^^Q, where n = N/L is the density of particles. For 
a mesoscopic system, the derivatives should be replaced by finite differences. In order to 
keep the magnetization fixed, the charge compressibility is then computed by adding both 
a spin up and a spin down particle to the ground state configuration; the explicit form of 
the derivative is 

N^ 

x;' = 2 X —[Eo{N + 2) + Eo{N - 2) - 2Eo{N)]. (4) 

The leading factor of two ensures that a pair of particles is introduced. 

Shastry and Sutherland^ have shown that the bulk spin susceptibility Xs in the Hubbard 
model is related to the spin stiffness by D^Xs = l/27r^n^ due to a remarkable property of the 
BA equations @. Thus in the thermodynamic limit we only need to compute either spin 
stiffness or susceptibility. However, for a mesoscopic system, this relation does not hold due 
to an avoided level crossing at = 27r.0 We remark that this situation is very much as for 
the Heisenberg-Ising model in the momentum vr sector.lli Following Ref. we then write 
for the spin susceptibility 

1 N^ 

x:' = ^ [EoiN, M = iV/2 - 1) - Eo{N, N/2)] . (5) 

Note that in the M = N/2 — 1 sector the ground state has the quantum numbers In = 
-{N + l)/2 + n for 1 < n < and Jj = -{N/2 + l)/2 + j for 1 < j < M. 
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The behavior of both Xc and Xs has been reported previously in the thermodynamic 
limit.llli2l Our mesoscopic results are qualitatively the same and differ at the most by 15% 
from their thermodynamic values. Thus we refrain from including the corresponding figures 
here. However, we emphasize that these small deviations will be quite important in the 
following sections. 



2. Boson representation of the repulsive Hubbard model 

Away from half filling, the low-energy and large distance behavior of a one-dimensional 
fermion system with spin-independent interactions is described by the Hamiltoniani'i 

H = H, + Hs + I dx cos iV8(j)six)), (6) 



where 



^Jdx [(t;,ir,)(7rn,)2 + (^-^) (7) 



for r = c,s. This Hamiltonian describes the most general ID Hamiltonian with spin con- 
serving interactions, provided that the proper values for and v^. are used.ill The c and 
s parts of the Hamiltonian describe the charge and spin degrees of freedom of the system 
respectively. The operator H^ is the momentum density conjugate to (pr and these operators 
obey Bose-like commutation relations: [Ilr{x), (f)r{x')] = —i5r,^5{x — x'). a is a short-range 
cutoff parameter of the order of the lattice constant. The gi_<^ term represents scattering 
between opposite spins with a momentum transfer close to 2k f. The umklapp scattering 
transferring two particles from —kf to kf, involves a momentum transfer Akf = 2tt which 
in the half filled band {kf = tt/2) corresponds to a reciprocal lattice vector. Away from 
the half-filled case this term does not contribute and therefore has not been included in the 
Hamiltonian (^. The case gi± = describes independent long- wavelength oscillations of 
the charge and spin densities with linear dispersion relation uJr{k) = Vr\k\. For gij_ ^ the 
cosine term has to be treated perturbatively. 
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For the Hubbard Hamiltonian the values for Kc and Vc are given by the following 
relations, 



Here Dc and Xc are defined as in Eq. (|I|) and (^), respectively, and n = N/L is the density 
of particles. For gij_ = 0, we can similarly write 



with Ds and Xs defined as in Eq. (0) and and the leading factor of 2 is due to our 
definition of the spin flux However, Dg and Xs have been calculated from the exact 
solution of the Hubbard model, e.g., including the 2k f scattering terms represented by the 
gi± term in the bosonized description of Eq. (^. Therefore, the identifications of Eq. 
are not exact and only valid perturbatively in gij_. Since for the Hubbard model gi± = U, 
this means that Eq. is valid only for small U. 

Let us briefly elaborate on this point: For a macroscopic system gi± renormalizes to 
zero.ii The spin stiffness in this limit is given as Dg = -^VsK* where K* = 1 is the 
renormalized value of Ks- However, for a finite system, one generates additional terms 
of 0{gi±) coming from the term proportional to gi±. We therefore define the spin stiffness 
as Dg = -^VgKs with Vg and Kg given as in Eq. (^). When Kg = 1, we find the simple 
relation between the spin stiffness and the spin susceptibility as in last section. However, 
for mesoscopic systems deviation of Kg from its macroscopic fixed point value occurs as a 
consequence of the avoided level crossing mentioned earlier. We find that Kg varies from its 
non interacting value Kg = 1 up to Kg 1.15 for f/ = 20 in a ring of size L = 100. 




(8b) 



(8a) 




(9b) 



(9a) 
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B. The attractive interaction regime 



The attractive Hubbard model is again described by the Hamiltonian of Eq. but 
now U < 0. The attraction gives rise to Cooper- hke bound pairs of spin up and spin down 
particles. These pairs scatter off each other without diffraction. Hence for even and 
zero magnetization, we can write down effective BA equations for M = N/2 pairs.i'ii Let 
a{ki) = 2arcsin (jvi + i^^. Then the BA equations for the pairs are given as, 

M 



L Re{a{ki)/2} = 2TcJi + 2$p + ^arctan 



77 ^i) 



(10) 



where $p = $-f = $| is the flux acting on the Cooper pair. The factor of two in front 
of <l>p arises due to the charge of the Cooper pair. The quantum numbers for the ground 
state are AJ^ = (M + l)/2 mod 1, and the corresponding ground state energy is given 
as Eo{^p) = — 4X]fii cos[Re{a(fci)/2}] cosh[Im{Q;(/ci)/2}]. There exists a gap in the spin 
excitation spectrum equal to the binding energy of a pair and only the charge sector remains 
gapless. Hence only the charge sector in the bosonized Hamiltonian is retained and the only 
relevant parameters in (|^) are and Vc- In order to make the reference to the Cooper pairs 
explicit, we shall rename these parameters Kp and Vp. 

The pair stiffness is defined as before Dp = {L/2) (9^i?o($p)/<9$p and the pair compress- 
ibility is computed by changing the ground state configuration by adding a Cooper pair. 
The thermodynamic expression for the pair compressibility is Xp ^ = {m? / L)d'^ Eq{M) / dm? , 
where m is the density of Cooper pairs. The explicit form of the derivative for a mesoscopic 
system is Xp^ = (M^/L) [Eo{M + 1) + Eo{M - 1) - 2Eo{M)]. We can now compute Kp 
and Vp from the ground state energies. The relations are 



Kp = Tim, 



'""^ (lla) 



\Xp' 



= —\ ' ■ lib 
m V 2 

The error involved in using the Cooper Pair representation is of 0(exp[— 2f/L]). From 
Fig. we see that down to \U\ = 0.01 there is very good agreement of the pair stiffness in 
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the attractive regime with the charge stiffness found in the repulsive regime. This justifies 
our use of the effective BA equations even for such small U. The general features of the pair 
stiffness are the same as that seen for the charge stiffness. One sees a reduction in the pair 
stiffness as the attraction is increased. Note that for small \U\ and density n, the value of 
Dc is always larger or equal to Dp. Here the charge currents in the repulsive case are larger 
than in the attractive case. For densities close to half- filling and U > 1, this behavior is no 
longer true as Akf scattering terms become dominant in the repulsive regime. 

IV. EFFECT OF DISORDER 

In the limit of weak disorder, the interaction between the particles and disorder can 
be parameterized by two uncorrelated Gaussian random fields t] and ^0 These two fields 
describe the forward and backward scattering by the impurities. The forward scattering 
term can be treated exactlyil^ in one dimension and is found not to contribute to the 
conductivity. However, the effect of backward scattering is very important and leads to 
localization in the non-interacting limit. ^ and ^* correspond to the part of the random 
potential which has Fourier components close to 2k f. Higher Fourier components are less 
effective and do not correspond to low energy processes. A notable exception is the 4k f term 
which we will have reason to discuss in the last section. In terms of the boson variables the 
impurity coupling to the particle density is given as 

H^ = Jdx e*^*^(^) cos (720,,) + h.c, (12) 

where ^ is a Gaussian with = A^6{x — x'). Contrary to the free case the charge 

and spin degrees of freedom are no longer independent but are coupled through the random 
potential. 

Generally both impurity backscattering and the interaction term gij_ give rise to divergent 
terms in a perturbation calculation. Hence, a perturbative approach in the disorder and 
in the interaction ^fi^ is used to generate the renormalization group equations under a change 
of the length scale a e'a,0@ where a is the lattice spacing. The equations are 
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dK^{l) _ 1 
dl ~ ^2 

dKsjl) _ 1 r 
dl 2 



A(/), (13a) 



\2 T^2 



A{l)+y{iy Ki, (13b) 



dvc{l) vlK, 



dl 2vs 



dl 2 
dyil) 



" ^A(/), (13c) 
A(/), (13d) 



dl 
dA{l) 

dl 



[2-2KS)]y{l)-m, (13e) 
[3-KS)-Ksil)-yil)]A{l), (13f) 



with the dimensionless quantities A and y defined as: A = {2A^a/nv^) [vs/vc]^" and y = 
Qix/tcVs- The development of the RG equations is hmited to first order in A^ and to second 
order in gi± and hence the renormahzations of Kr and Vr can be neglected on the right hand 
side of the first three equations. 

For a mesoscopic system of size L, the infra-red cutoff is expected to be given by L s.t. 
L ~ e'a, or equivalently, / ~ InL/a. Thus we may calculate charge and spin stiffness at 



finite size by using the values of Kc{l), Ks{l),Uc{l) and Cs{l) and the formulas of section |T1 



The RG equations for finite disorder have to be numerically integrated taking appropriate 
initial conditions at / = 0. Using perturbative values as initial data for the RG equations, 
e.g., Vs = Vf and Ks = 1 + y /2, we see that we get an increasing value for the spin stiffness 
Ds = -^VsKs with increasing repulsion. This does not agree with the exact Bethe Ansatz 
solution for the spin stiffness as shown in Fig. |^ which clearly shows a reducing stiffness. 
Therefore, we shall use the clean Hubbard model parameters that were computed in section 
pni and gi±_ = U . In the following sections we discuss the results of the integration primarily 
for n = 0.3. The results for fillings corresponding to M = 10, 22, 42, 50, 62, 70, 82, 90 and 96 
are qualitatively the same. 
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A. The repulsive case 



For the repulsive case we find that A always fiows to oo as shown in Fig. ^. Assuming 
that there is no other fixed point at intermediate coupling, this whole region can then be 
identified with the localized phase. The magnetic properties of the system will then 
depend of the renormalized value of y: 

(i) Increasing the repulsive interaction away from U = 0, we find that y fiows to — oo as 
shown in Fig. ^. The physical state corresponds to a non-magnetic system of localized pairs 
of spins which is equivalent to a charge density wave pinned by the impurities (PCDW).ilii 
In Fig. H we see that the disorder scales less rapidly to infinity as the interaction strength is 
increased. This implies that pinning of the charge density wave is harder in the presence of 
a repulsive interaction. Thus we expect the interacting charge currents to be enhanced by 
the disorder. However, as shown in section |II|, there is an initial reduction in the interacting 
charge currents to due the presence of the lattice. In Fig. ^ we show the renormalization 
of the charge stiffness for a fixed filling n = 0.3 and disorder strength for various values 
of the interaction strength. We see that there is a crossover between the non interacting 
and interacting stiffnesses for finite 1^, or, equivalently, finite disorder A^(/c, U) and we find 
enhancement only for disorder values larger than A^{lc, U). The crossover regions for a fixed 
value of the disorder as a function of n and U are shown in Fig. |^. Previous studies of 
the disordered Hubbard rin g§S treat the Hubbard interaction in first order perturbation 
theory for finite disorder and thus cannot identity this crossover. 

In Fig. H we further note that y scales less rapidly to — oo as the interaction strength 
is increased. As noted earlier, the physical state corresponds to a system of localized spins 
and hence larger the U value lesser the localization. We therefore find that the spin currents 
are less pinned by the impurity as the repulsion is enhanced. In Fig. ^ we again find due to 
the reduced pinning effect there exists a critical disorder strength A^{ls,U) where there is 
a crossover between the non-interacting and interacting spin currents. Ig here is the length 
at which the crossover occurs and A^{ls, U) is the value of the disorder at Is for a given U . 
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Note that we always find Ic < h {^^{h, U) < A^{ls, U)). We remark that the actual value 
of A^{ls,U) is independent of particle density n. However, we expect a strong dependence 
on the magnetization. As we have restricted our study to the sector of zero magnetization, 
we do not observe this later dependence here. 

(ii) In Fig. ^ we see that further increase of U results in y flowing to +00. Hence there 
is a strong repulsion of up and down spins, and the particles start to localize as isolated 
spins on randomly distributed sites. The magnetic properties have been identified earlier as 
being typical of a random antiferromagnet (RAF). 00 

The charge currents in the RAF phase will continue to be enhanced. In Fig. ^, we see 
that, e.g., the charge stiffness for n = 0.3 and U values of 0.8 and 1.0 is still above the 
non-interacting stiffness. Fig. ^ clearly shows that these values of U already belong to the 
RAF fixed point y* = 00. 

The strong repulsion of up and down spins, however, gives rise to a drastic fall in the spin 
currents as seen again for n = 0.3 in Fig. |^. We thus observe no enhancements in the spin 
currents in the RAF phase. This may be used to distinguish the two different phases of the 
system for repulsive U. In Fig. |^ the region above the line represents the RAF state where 
no enhancement of the spin currents is found, whereas the region below the line corresponds 
to the PCDW state. The reader should compare Fig. | with Fig. ^ The maxima of each 
curve in Fig. ^ corresponds to the transition point between PCDW and RAF phases as in 
Fig. p. We note that for values of n and U such that the fixed point of the system belongs 
to the RAF phase, already small disorder values will localize the spins and thus lead to an 
enhancement of the charge stiffness. 



B. The attractive case 



For negative U there exists a gap in the spin excitation spectrum and only the charge 
sector remains gapless. The RG Eqs. (0) can be reduced by taking into account only the 
charge (i.e. pair) excitations,0'ii 
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A(Z), (14b) 



^ = -^ir|A(/), (14a) 

dvpjl) _ v^Kp 
dl ~ ~ 2 

^ = [3-Kp{l)]A{l), (14c) 
with A(/) = {2CsA^a/TTVp) and Cs a constant of order unity. These equations have been 



studied previously and we find in agreement with Ref. |T0|: From Eq. ( |14c| ) we see that as 
Kp < 3, the disorder will always scale to large values thereby always leading to localization 
of the particles. This implies for the second equation that Vp —>■ 0. The stiffness therefore 
always reduces in the presence of disorder. As the attraction is increased the disorder 
scales faster to infinity and since in the attractive Hubbard model the ground state contains 
strong charge density fiuctuations they get easily pinned to the disorder. Due to the large 
reduction in the bare stiffness and an increased pinning effect to the disorder we therefore 
see no enhancement in the persistent currents as the attraction is increased. 



V. CONCLUSIONS 

We have studied in this work the behavior of charge and spin stiffness constants of a ID 
disordered mesoscopic Hubbard ring. For the attractive Hubbard model without disorder, 
we find in agreement with previous studies0 that the charge stiffness gets reduced with 
increasing \U\ and is always smaller than the corresponding charge stiffness for the repulsive 
regime at densities not close to half-filling. For finite disorder, we too find no enhancements 
in the stiffness for finite \U\ over the non interacting currents. 

For the repulsive Hubbard model without disorder, we again see that the charge and spin 
stiffnesses get reduced as U is increased. Additionally, we observe a very strong reduction as 
we approach the half-filled situation for finite U. However, for finite disorder, the situation 
is very different from the attractive case: For small U, the inclusion of disorder drives 
the system into a localized state with enhancement of both charge and spin stiffnesses. 
However, for a mesoscopic system, enhancement is only observed for large enough values of 
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the disorder, or, equivalently, large enough system sizes. The physical state in this phase 
corresponds to a non-magnetic system of localized pairs of spins which is equivalent to a 
charge density wave pinned by the impurities (PCDW).^^ As the interaction is increased 
the tendency towards pinning is found to reduce as the disorder scales less rapidly towards 
oo. 

For larger U, the effective interaction between up and down spins becomes repulsive 
and we identify the physical state as corresponding to isolated spins localized on randomly 
distributed sites (RAF). 00 This strong repulsion between unlike spins reduces the spin 
current drastically and we find no enhancements in the spin current. However, the charge 
currents in the RAF phase are still enhanced as compared to the non interacting current. 

We emphasize that the transition from PCDW state to RAF state at intermediate values 
of the interaction strength U can be identified by studying the different behavior of Dg in 
the two phases. 

In this study, we have neglected Umklapp and Akj scattering terms both in the bosonized 
formulation of the Hubbard Hamiltonian and in the RG treatment of the disorder. Thus 
processes close to half-filling and U ^ 1 are presumably not included in our study. In order 
to at least qualitatively see what happens for large U, we have studied a 4 site Hubbard ring 
by exact diagonalization. Then due to the finite value of Ic, we need to go to large values of 
the disorder strength > 5 in order to observe the crossover in D^. However, we still do 
find a crossover even for t/ = 100. 
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FIGURES 

FIG. 1. The charge Dc and pair Dp stiffness for the Hubbard ring of size 

L = 100, as a function of the interaction strength for different values of the fiUing 
N = 2M = 10, 22, 30, 42, 50, 62, 70, 82, 90, 96 from bottom to top. 

FIG. 2. The spin stiffness Dg for the Hubbard ring of size L = 100, as a function of the 
interaction strength for different values of the fUhng N = 2M = 10, 22, 30, 42, 50, 62, 70, 82, 90, 96 
from bottom to top. Note that for large and U, it becomes increasingly difficult to distinguish 
different Dg curves due to limited numerical resolution. 

FIG. 3. Values of the disorder A as a function of the system size / ~ InL/a 



obtained by numerically integrating Eq. |13| for different values of the interaction strength 
U = 0,0.01,0.05,0.1,0.15,0.2,0.3,0.4,0.5,0.6,0.8, 1.0, 1.2,1.4, 1.5, 1.6, 1.8, 2.0, 4., 6., 8., 10., 15 and 
20 (small U is represented by small dashes) for a fixed density n = 0.3. The initial value for the 
disorder is fixed to A/vf = 5 x 10~^. The disorder scales less rapidly to infinity as the interaction 
strength is increased. 

FIG. 4. Values of y as a function of the system size / ^ InL/a. Density and initial disorder 
strength are as in Fig. ^, U varies from to 20. Increasing U is indicated by increasing dash 
length, y scales less rapidly to — oo as the interaction strength is increased from to 0.6. Starting 
at U = 0.8, y then scales to +oo. 

FIG. 5. Values of the charge stiffness Dc as a function of the system size / ~ InL/a. Density 
and initial disorder strength are as in Fig. ^ For clarity, U is only shown for values equal to 
0.0,0.1,0.2,0.4,0.6,0.8 and 1. Increasing U is indicated by increasing dash length. There exists 
a crossover region beyond which enhancement of the currents is observed. Note that the lines for 
U = 0.8 and U = 1 correspond to the RAF phase and are well above the U = value. 
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FIG. 6. We plot the disorder strength A^(/c) = {'n'v1^/2a)[vc/vs\^'' at which the enhancement 
of interacting over non interacting Dc occurs as function of interaction strength U . Different curves 
correspond to different fillings, i.e., N = 2M = 10,22,30,42,50,62,70,82,90,96 from bottom to 
top. Note that we have only limited resolution in U corresponding to the possible values of the 
interaction strength as in Fig. ^. 

FIG. 7. Values of the spin stiffness Ds as a function of the system size / ~ InL/a. Density 
and initial disorder strength are as in Fig. ^, U varies from to 20. Increasing U is indicated by 
increasing dash length. There exists a crossover region beyond which enhancement of the currents 
is observed for < C/ < 0.6. After the PCDW-RAF transition between Q.Q <U < 0.8, Dg rapidly 
scales to zero. 

FIG. 8. Using the qualitatively different behavior of Dg, we identity the boundary between 
the PCDW and RAF phases. The data points represent the largest values of U for a given filling 
at which we still see a crossover in Dg and are thus part of the PCDW phase. Note that we have 
only a limited resolution in U corresponding to the possible values of the interaction strength as 
in Fig. ^. 
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